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ABSTRACT 

We present measurements of the dimensionless spin parameters and inner-disk inclination of 
two stellar mass black holes. The spin parameter of SWIFT J1753.5-0127 and GRO J1655-40 
are estimated by modelling the strong reflection signatures present in their XMM-Newton 
observations. Using a newly developed, self-consistent reflection model which includes the 
blackbody radiation of the disk as well as the effect of Comptonisation, blurred with a rel- 
ativistic Une function, we infer the spin parameter of SWIFT J1753.5-0127 to be Q.lQtlW- 
The inclination of this system is estimated at SSly degrees. For GRO J1655-40 we find that 
the disk is significantly misaligned to the orbital plane, with an innermost inclination of 301^q 
degrees. Allowing the inclination to be a free parameter we find a lower limit for the spin of 
0.90, this value increases to that of a maximal rotating black hole when the inclination is set 
to that of the orbital plane of J1655-40. Our technique is independent of the black hole mass 
and distance, uncertainties in which are among the main contributors to the spin uncertainty 
in previous works. 
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1 INTRODUCTION 

Black holes can be described by two observable parameters, mass 
and spin. To date there are over twenty stellar mass black holes 
with dynamically constrained mass (for a review see McClintock 
& Remillard 2006); however for just a handfuU of these systems do 
we have any measurements of their spin. 

Shafee et al. (2006) and McClintock et al. (2006) have re- 
ported values of the dimensionless spin parameter, a for 4U 1543- 
47 of 0.7-0.85, GRO J1655-40 of 0.65-0.75, and GRS 1915+105 
of 0.98-1. More recently Liu et al. (2008) reported a value of 
0.77 ± 0.05 for M33 X-7. Their approach relies on modelling the 
thermal continuum seen in the thermal X-ray spectra and requires 
sources to be selected in the high-soft state (for a review on spectral 
states see McClintock & Remillard 2006. Any powerlaw emission 
is then minimal and the spectrum resembles a quasi-blackbody con- 
tinuum. Furthermore, precise measurements of the mass and dis- 
tance to the black hole, as well as the inclination of the system, are 
essential. 

We have recently reported precise measurements of the spin 
parameter in GX 339-4 (Reis et al. 2008; Miller et al. 2008a) using 
the reflection signatures in the spectrum. These reflection features 
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arise due to reprocessing of hard X-ray by the cooler accretion disk 
(Ross & Fabian 1993), and consist of fluorescent and recombina- 
tion emission lines as well as absorption features. In the inner re- 
gions of an accretion disk the various reflection features become 
highly distorted due to relativistic effects and Doppler shifts. The 
shape of the prominent ¥e-Ka fluorescent line, and more impor- 
tantly the extent of its red wing, can give a direct indication of the 
radius of the reflecting material from the black hole (Fabian et al. 
1989, 2000; Laor 1991). The stable circular orbit around a black 
hole extends down to the radius of marginal stability, rms, which 
depends on the spin parameter (e.g Bardeen et al. 1972). A major 
advantage of using reflection features to obtain the spin of the black 
hole is that these features are completely independent of black hole 
mass and distance and can thus be used for any system where both 
or either of these parameters are unknown (for a review see Miller 
2007). For GX 339-4 we found that for both the very-high and the 
low-hard state the accretion disk extends to the innermost stable cir- 
cular orbit, Tms at a radius of ~ 2.05rg, where rg= GM/(? (Reis 
et al. 2008; Miller et al. 2008a). This implies a spin of ^ 0.935 for 
GX 339-4. 

In this paper we use the method adopted by Reis et al. 
(2008) which uses a specially developed spectral grid, REFHI- 
DEN (Ross & Fabian 2007), to obtain the spin parameter of SWIFT 
11753.5-0127 and GRO 11655-40. 11753.5-0127 was first detected 
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in hard X-rays by the Burst Alert Telescope (BAT) on the SWIFT 
satellite on 2005 May 30 (Palmer et al. 2005). Using RXTE and 
XMM-Newton data, Miller, Homan & Miniutti (2006, hereafter 
M06) showed the presence of a cool (kT ~ 0.2 keV) accretion disk 
extending close to the rms in the low-hard state (LHS) of the sys- 
tem. The presence of this cool accretion disk was later confirmed 
by Ramadevi & Seetha (2007) and more recently by Soleri et al. 
(2008) using multiwavelength observations of the source. 

GROJ1655-40 was discovered by the Burst and Transient 
Source Experiment (BATSE) on-board of the Compton Gamma 
Ray Observatory (CGRO) on 1994 July 27 (Zhang et al. 1994). The 
mass of the compact object has been dynamically constrained to 
> 6.0 Mq (Orosz & Bailyn 1997). ASCA observation of J1655-40 
from 1994 through 1996 provided the first clear detection of 
absorption-line features in the source (Ueda et al. 1998; see also 
Miller et al. 2008b). Using archival ASCA data, Miller et al. (2004) 
showed evidence of highly skewed, relativistic lines, and suggested 
an inner-radius of reflection of ~ l-4rg which would indicate a 
rapidly spinning black hole. This was later confirmed by Diaz- 
Trigo et al. (2007, hereafter DT07) using simultaneous XMM- 
Newton and INTEGRAL observations of J 1655-40 during the 2005 
outburst. The high spin suggested by these authors is in contrast 
with the spin parameter reported by Shafee et al. (2006) of 0.65- 
0.75 using the thermal thermal continuum method. In the following 
section we detail the observation, analyses procedure and results. 



2 OBSERVATION AND DATA REDUCTION 

5iy/Fr J1753.5-0127 was observed in its low-hard state by XMM- 
Newton for 42 ks, starting on 2006 March 24 16:00:31 UT and si- 
multaneously by RXTE for 2.3 ks starting at 17:26:06 UT (M06). 
The EPIC-pn camera (Struder et al. 2001) was operated in "timing" 
mode with a "medium" optical blocking filter. For GRO J 1655-40, 
observations were made by XMM-Newton for 23.9 ks on 2005 
March 18 15:47:13 (hereafter Obs 1) and again on 2005 March 
27 08:43:59 (hereafter Obs 2) for 22.3 ks (DT07). The source was 
found to be in the high-soft state and was observed with the EPIC- 
pn camera in "burst" mode with a "thin" optical blocking filter. 
Starting with the unscreened level 1 data files for all aforemen- 
tioned observation we extracted spectral data using the latest XMM- 
Newton Science analyses System v 7.1.0 (SAS). For J1753.5-0127 
events were extracted in a stripe in RAWX (20-56) and the full 
RAWY range. RAWX 30-43 and RAWY 5-160 was used for 
J 1655-40. For both sources, bad pixels and events too close to chip 
edges were ignored by requiring "FLAG = 0" and "PATTERN^ 4". 
The energy channels were grouped by a factor of five to create 
a spectrum. Background spectra were extracted for both sources 
from an adjacent region of similar RAWY and RAWX. In the case 
of J1655-40 the background is negligible due to the high source flux 
and it was not used for the analyses that follows. In both cases, indi- 
vidual response files were created using the SAS tools RMFGEN and 
ARFGEN. The EPIC-pn resulted in a total good-exposure time of 
40.1 ks for J1753.5-0127. Due to the low duty-cycle of the "burst" 
mode (3%) the total good exposure time for J1655-40 is ~ 0.7 and 
0.6 ks for Obs 1 and 2 respectively. RXTE data for J1753.5-0127 
were reduced in the standard way using the HEASOFT v 6.0 soft- 
ware package. We used the "Standard 2 mode" data from PCU-2 
only. Standard screening resulted in a net Proportional Counter Ar- 
ray (PCA) and High-Energy X-Ray Timing Experiment exposures 
of 2.3 and 0.8 ks respectively. To account for residual uncertainties 
in the calibration of PCU-2, we added 0.6 per cent systematic error 



to all energy channels. The HEXTE-B cluster was operated in the 
"standard archive mode". The FTOOL GRPPHA was used to give at 
least 20 counts per spectral bin. 

We restrict our spectral analyses of the XMM-Newton EPIC- 
pn data for J1753.5-0127 to the 0.7-10.0 keV band. For J1655-40 
we use 0.7-9.0 keV due to the uncertain calibration above 9.0 keV 
for the EPIC-pn burst mode. The PCU-2 and HEXTE spectrum 
are restricted to 5.0-25.0 and 20.0-100.0 keV band respectively. A 
Gaussian line at 2.2-2.3 keV is introduced when fitting the EPIC- 
pn spectrum due to the presence of an instrumental feature in this 
energy range that resembles an emission line. This feature is likely 
to be caused by Au M-shell edges and Si features in the detector. 
All parameters in fits involving different instruments were tied and 
a normalisation constant was introduced. XSPEC v 12.4.0 (Arnaud 
1996) was used to analyse all spectra. The quoted errors on the 
derived model parameters correspond to a 90 per cent confidence 
level for one parameter of interest (Ax^ = 2.71 criterion). 



3 ANALYSIS AND RESULTS 

Fig. 1 shows the data/model ratio for J1753. 5-0127 fitted with a 
simple absorbed powerlaw (PHABS model in XSPEC with A^h 
fixed at 1.72 x 10^^ cm"^) in the energy range 2.5-5.0 and 8.0- 
10.0 keV and then extended to the full energy range. As first noted 
by Miller, Homan & Miniutti (2006) and more recently confirmed 
by Soleri et al. (2008), J1753. 5-0127 shows clear indications of the 
presence of a cool accretion disk with a temperature of ~ 0.2 keV 
in its low-hard state. The presence of an excess at ~ 6.9 keV is 
indicative of Fe-A'a emission. Similarly for J1655-40, we fit the 
0.7-4.0 and 7.0-9.0 keV energy range with a power law and a fur- 
ther multi-color disk blackbody (MCD; Mitsuda et al. 1984) modi- 
fied by absorption in the interstellar medium. The value of A'h was 
kept constant for the two observations. The data/model ratio for 
J 1655-40 extended to the full energy range is shown in Fig. 2. The 
presence of a broad Ve,-Ka line emission extending to just over 
3.0 keV as well as absorption features are clearly seen. Fitting the 
low/hard state of J1753. 5-0127 with a simple absorbed powerlaw 
results in an unacceptable fit (Table 1) with various residuals in 
both the soft and the ¥e,-Ka energy range. 

The presence of a soft disk excess and broad ¥e,-Ka emission 
line is usually modelled phenomenologically with a combination 
of a multi-color disk component such as DISKBB (Mitsuda et al. 
1984) and a relativistic line such as the LAOR line profile (Laor 
1991). This combination of components can give robust spin mea- 
surements when the presence of the ¥s-Ka emission line is signif- 
icantly above the continuum. However it should be noted that it is 
only an approximation since the identification of rin, as determined 
from LAOR assumes a hard wired spin parameter of a = 0.998. 
Furthermore, the LAOR line profile is a phenomenological model 
for a relativistic emission line. The effects that extreme gravity have 
on the reflection signatures are not limited to the Fe-A'a line profile 
and thus a more thorough method to constrain the spin would in- 
volve modelling all of the reflection signatures present in the spec- 
tra. In what follows, we use the reflection model developed by Ross 
& Fabian (2007, REFHIDEN) to model all the reflection signa- 
tures as well as the disk blackbody emission in a self-consistent 
manner. This approach was detailed in Reis et al. (2008) for the 
galactic black hole GX 339-4. 

The reflection features produced by the illuminated surface 
layer of an accretion disk are largely dependent on the ionisation 
state of the disk, and thus an important quantity is the ionisation 
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Figure 1. Data/model ratio for J1753. 5-0127 obtained by fitting the en- 
ergy range 2.5-5.0 and 8.0-10.0 kcV with an absorbed powerlaw. It is clear 
that a semi-blackbody component as well as an iron reflection signature is 
present in the low-hard state of J1753. 5-0127. The data have been rebinned 
for plotting purposes only. 
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Figure 2. Data/model ratio for J1655-40 using a simple absorbed power law 
and disk blackbody fitted in the energy range 0.7^.0 and 7.0-9.0 keV. Obs 
1 and 2 are shown in black and red respectively. The presence of a broad 
¥e.-Ka line and complex absorption is clearly seen. The data have been 
rebinned for plotting puiposes only. 



parameter ^ = AnFh/ Hden, where Fh is the hard X-ray flux il- 
luminating a disk with a hydrogen density Hden (Matt, Fabian & 
Ross 1993). The REFHIDEN reflection model incorporates the 
importance of the thermal emission from the disk mid-plane in de- 
termining the ionisation state and thus reflection features. Whereas 
previous reflection models such as the Constant Density Ionised 
Disk (CDID; Ballantyne, Ross & Fabian 2001) and REFLIONX 
(Ross & Fabian 2005) vary over the ionisation parameter ^ and does 
not account for thermal emission, in REFHIDEN both the number 
density of hydrogen, //den, and the ratio of the the total flux il- 
luminating the disk to the total blackbody flux emitted by the disk 
(-Riiiujn/flfl) are implicit parameters. The value of kTfox the black- 
body entering the surface layer from below and the power-law pho- 
ton index are further parameters in the model. The disk reflection 
spectra is convolved with relativistic blurring kernel KDBLUR, 
which is derived from the code by Laor (1991). The parameters for 
the blurring kernel are the emissivity index gi„ , where the emissiv- 
ity (er) is described by a powerlaw-profile such that e,. = j-^^in^ 
disk inclination i, the inner disk radius nn, and the outer disk radius 



Table 1. Results of fits to XMM-Newton EPIC-pn data for J1753.5-0127 
with both a simple absorbed powerlaw and the self-consistent reflection 
model REFHIDEN. 



Parameter 



Simple 



REFHIDEN 



Nn (102 



r 

kTikeV) 

i/donCxlO^^Hcm- 

Rlllum/BB 

A^REFHIDEN 

gin 

nn(rg) 
(deg) 



1.54 ±0.13 
1.631 ±0.003 
0.0591 ± 0.0002 



1.52 
1.54 
0.050 
0.193 
1.40 



2191.9/1859 



0.02 
0.03 
+0.01 
0.04 
+0.003 
0.004 
+0.003 
0.002 
+0.13 
0.27 
Q fi+0-4 
y-0_0.6 

+0.0004 
0.0006 
+6 
1 

+0.7 
0.6 
+2 
7 



0.0013 
4 
3.1 

55 

1839.5/1852 



Notes.- The self-consistent model is described in XSPEC as 
PHABSx KDBLUR x(PL+REFHIDEN). The normalisation of each 
component is referred to as A'. All errors refer to the 90% confidence range 
for a single parameter. 



which was fixed at 400rg. The power law index of REFHIDEN is 
the same as that of the hard component. 



3.1 Self-consistent reflection and disk emission 

• SWIFT J1753.5-0127: The model provides an excellent fit to 
the data for J1753.5-0127 with x^/f = 1839.5/1852. Table 1 de- 
tails the parameter values with all errors corresponding to the 90 
per cent confidence range. The data/model spectrum is shown in 
Fig. 3a with the fit extended to 100.0 keV using RXTE data shown 
in the inset. A normalisation constant was added to account for 
flux mismatch between the instruments. The best fit model prior 
to gravitational blurring is shown in Fig. 3b. The value of the inner 
radius obtained from the gravitational blurring of the reflection fea- 
tures is constrained to be rin= S.l^Q grg. The strong constraint on 
the value of nn can be better appreciated in Fig. 4a, where the 90 
per cent confidence level is shown in the plot obtained with the 
"steppar" command in XSPEC. Fig. 4b shows a similar constraint 
obtained for the inner accretion disk inclination of J1753. 5-0127. 

As can be seen from Table 1, the value of the emissivity in- 
dex, Qin has been poorly constrained. A value of 3 is expected for 
a standard accretion disk (Reynolds & Nowak 2003) with steeper 
values usually being interpreted as emission from a small, compact, 
centrally located X-ray source, such as expected from the base of 
a jet (Miniutti & Fabian 2004). In order to investigate any degen- 
eracy between the value of the inner radius and the unconstrained 
emissivity we explored their parameter space using the "contour" 
command in XSPEC. All parameters except for F and Riuum/BB 
were free to vary. Fig. 5a shows the 68, 90 and 95 per cent con- 
tour for two parameters of interest. Although the emissivity index 
is poorly constrained, it can be seem from Fig. 5a that the value 
of the inner radius is not strongly affected by this uncertainty. In 
what follows we will thus freeze the value of the emissivity at 
the best fit value shown in Table 1. Fig. 5b and 5c shows simi- 
lar contour plots for both inclination and Hden versus inner radius 
respectively. For a large range of inclination and Hden the value 
of r,n remains approximately between 2^rg with a best fit value 
of approximately 3rg, in accordance to that shown in Fig. 4a. As- 



suming that nn= 3.1 



(Fig. 4a) is the same as the radius of 
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Figure 3. (a): Data/model ratio for J1753. 5-0127 in its low-hard state. The model assumes a powerlaw emissivity profile and constitutes of a powerlaw and 
the disk reflection model REFHIDEN. The inset shows the model extended to 100.0 keV using RXTE PCA (red) and HEXTE (green). Data have been 
rebinned for plotting purposes only, (b): Best-fit model prior to gravitational blurring showing the reflection features. The total model, powerlaw and reflection 
components are shown in black, blue and red respectively. 




Figure 4. (a): r-jn plot for J1753. 5-0127. A value of rin= S.l^Q grg is found at the 90 per cent confidence level for one parameter of interest 

(Ax^ = 2.71 ciiterion) shown by the solid hoiizontal line, (b): Similar plot for the disk inclination. A value for the disk inclination of 55^^ is found at the 
90 per cent confidence level for one parameter of interest. 
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(c) 

Figure 5. (a): Emissivity versus inner radius contour plot for J1753. 5-0127. The 68, 90 and 95 per cent confidence range for two parameters of interest are 
shown in black, red and green respectively, (b): Similar plot for the inclination versus inner radius and (c): Hydrogen density versus inner radius. It can be 
seen that for the full range of the emissivity, inchnation and hydrogen density, the inner radius is constrained between approximately 2-Arg at the 90 per cent 
confidence level for two parameters. 
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Figure 6. Plot of the innermost stable circular orbit vs dimensionless spin 
parameter (solid curve). The constraints imposed by the innermost radius 
obtained for J1753. 5-0127 are shown by the intersection of the solid Une 
with the vertical region. The solid region shows similar constrains for the 
lower limit of the spin parameter in J1655-40. Based on this analysis, a 
dimensionless spin parameters of 0.76^Q'i[5 is found for J1753. 5-0127 and 
a lower limit of « 0.90 is found for J1655-40. 

marginal stability r-ms, we obtain a dimensionless spin parameter 
of 0.761 o ls for J1753.5-0127 in its low-hard state (Fig. 6). 
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Figure 7. Top: Best-fit data/model ratio for J1655-40 using Model 1 (Ta- 
ble 2). The various absorption features are clearly seen. Bottom: Similar 
as above for Model 2 covering the energy range 0.7-6.4 kcV. The model 
clearly results in a good fit for the red-wing of the \mn-Ka fluorescent 
line. Obs 1 and 2 are shown in black and red respectively. The data have 
been rebinned for plotting purposes only. 



• GRO J1655-40: Strong absorption features are clearly 
present in the spectra of J 1655-40 (see Fig. 2), and thus a fit with 
REFHIDEN should not immediately give an acceptable result. 
Fig. 7 (Top) shows the best-fit data/model ratio using REFHI- 
DEN. The various parameters for this model are described in Table 
2 (Model 1). It can be seen that, whereas in Fig. 2 there is evi- 
dence of both a broad line and various absorption features, this 
time the presence of the broad line has diminished (Fig. 7 Top). 
This fit constrains the inner radius to nn= l-SSlg q2 'g ^i^h 
/v = 4967.3/3296. The best fit using Model 1 seems to im- 
ply a relatively low inclination of 50°±1°, considerably less than 
the value of the binary inclination of 70.2°±1.9° (Greene, Bai- 
lyn & Orosz 2001). However, this value is still in agreement with 
that presented by Diaz-Trigo et al. (2007, see their Table I). The 
low inclination seen here could be due to the various absorption 
features masking themselves as the blue wing of an uon-Ka line 
profil4B 

Although various features contribute to the determination 
of the inner radius of emission and thus spin of the black hole, 
the most important of those features for the purpose of this work 
is the extent of the red wing in the ¥e.-Ka fluorescent line. In 
order to explore whether the absorption features are affecting the 
value of the inner radius, we froze the inclination at 70 degrees and 
modelled the spectra below 6.4 keV. The results of this fit is de- 
tailed in Model 2 (Table 2) and shown in Figure 7 (Bottom). By 
restricting the inclination to the known value of the binary incli- 
nation (~ 70 degrees) and fitting the spectra below 6.4 keV the 
innermost radius approaches that of a maximally rotating black 
hole with rin= 1.31 ± O.Olrg and xV^^ = 3230.0/2264. The 
majority of the contribution to is now coming from residuals 
between 0.9-1.1 keV possibly due to the photo-ionisation edge of 

^ The blue wing of the iron-Ka line profile can be used to determine the 
inclination of the source (see e.g. Reynolds & Nowak 2003; Fabian & Mini- 
utti 2005) 



O VIII at ~ 0.9 keV. In order to extend this fit to the full range 
we modelled the various absorption features in a phenomenologi- 
cal manner using two negative Gaussian absorptions (GABS model 
in XSPEC) fixed at energies of 6.7 and 6.97 keV corresponding to 
absorption of Helium-like Fe-XXV and Hydrogen-like Fe-XXVI 
respectively as well as absorption edges fixed at 8.8 and 9.3 keV. 
The optical depths r of the absorptions and edges were linked for 
self-consistency. The various parameters for this model covering 
the full energy range are detailed in Table 2 (Model 3). The value 
of the inner radius remains low (< 1.4rg) similarly to that of Model 
2. Allowing the inclination to vary over its full range improved 
the statistics with Ax'^ ~ —55.2 for one less degree of freedom 
(Model 4) and resulted in an inner radius of nn= 2.17lo i7rg 
(Fig. 8a) and inclination of 30°tio (^^8- ^^^^ spectra 

for JI655-40 including the absorption lines and edges is shown in 
Figure 9. Note that as we remove the restriction on the inclination, 
the emissivity profile of the two observations become less steep, 
with the indices constrained to qin = 3.87 ± 0.08 and 2.75^0 02 
for obs I and 2 respectively. Adding a possible nickel emission line 
at 7.48 keV with equivalent width of 30 and 17 eV for obs 1 and 2 
respectively resulted in an improved fit with Ax^ = —72.2 for 4 
extra degrees of freedom. As was the case for Model 2, most of the 
contribution to chi-squared are now coming from residuals in the 
soft energy range. By restricting this fit to above 1.5 keV we obtain 
X'/v = 3496.0/2957. 

Similarly to J1753. 5-0127 we investigate any degeneracy be- 
tween the inclination and inner radius in Model 4 by using the 
"contour" command in XSPEC. All parameters except for F and 
RiUum/BB were free to vary. Fig. 10 shows the 68, 90 and 95 
per cent contour for the two parameters of interest. It can be seem 
from Fig. 10 that for a large range of inclination the value of nn re- 
mains below 2.35rg with a best fit value of approximately 2.15rg, 
in agreement with the result presented in Fig. 8a. It is clear that the 
unknown inclination of the system plays an important role in the 
determination of the inner radius (and thus spin) of J1655-40, with 
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Table 2. Results of fits to XMM-Newton EPIC-pn data for J1655-40 with self-consistent reflection model REFHIDEN. 



Parameter 

r 

kT ( keV) 

Rlllum/BB 

N , 

Sin 

nn(rg) 
(' 

£^Gai,si(keV) 
a-( keV) 

T 

EEdgel (keV) 
£^Gabs2 (keV) 
cr(keV) 
r 

EEdge2 ( keV) 



Model 1 

Obs 1 

76+0 °l 
U- 'D_o.02 

2 59+0 03 
10 S+0-2 

Q 97+0.03 
-0.56 
n QQ+0.02 
'J-y'5_o.o4 
n Q+0.4 
0.0_o 2 



-0.3 



10.0. 

50 ± 1 



Obs 2 



O+0.03 
'-0.02 
1+0.3 

-0.1 



10.4: 

0.84 ±0.01 
10^ 

0.20:! 

0.29^ 



1+0.0 
-0.2 
1+0.03 
-0.01 
3+0.06 
-0.09 

4.5 ±0.1 



Model 2 
Obs 1 
0.759 ±0.001 
2.72 ± 0.08 

OQ + O.l 

y.u_o 2 
0.778 ± 0.009 
9.0l 
0.39^ 
1.4l 
10.0'' 



1+0.3 
-0.7 
1+0.05 
-0.10 
^+0.1 
-0.3 
1+0.0 
-0.3 

1.31 ±0.01 



70(f) 



Obs 2 



2 91+0.04 
^■^^-0.053 

11 ±3 
0.848lHj;^2 

0.6 ±0.2 
0.7 ±0.2 



10.0 



+0.0 



Model 3 
Obs 1 
0.765 ±0.001 



2.54' 



-0.02 
-0.01 



10.4^ 



^+0.1 
-0.2 

0.776+«-OO2 



7.4" 



■0.008 
■0.3 



4967.3/3296 



3230.0/2264 



-0.1 

11+0-01 

'-'•1^-0.03 

0.46t«:g^ 
10.0^0;° 
1.38 ±0.01 
70(f) 
6.7(f) 
0.08 ±0.02 
030+0 003 

U.UJU_o.o05 
8.8(f) 
6.97(f) 
0.06 ±0.02 

041+0 004 
U.U41_0 005 

9.3(f) 
4425.3/3289 



Obs 2 



2 99+0.01 

9 46+0-23 
^■™-0.02 
n yon+O.OOO 
0.78U_o 001 

lot; 

0.01+0-20 
0.02+0-10 

7 7c:+0.06 
' . ' on 



6.7(f) 
0.003 ±0.001 
0.20 ±0.06 
8.8(f) 
6.97(f) 
< 0.8 
n 014+0.016 

U.Ui^_0 008 

9.3(f) 



Model 4 
Obs 1 

„ 742+0-004 
'^^-0.002 
2 54+0-06 

o 1O+0.16 
».i»_0.12 

0.804 ± 0.004 

1+0.0 



10.010.6 

1+0.03 
-0.02 
,+0.01 
-0.04 



i.o: 

3.42:! 

3.87 ± 0.08 
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^■-^'-0.17 

6.7(f) 
0.09 ± 0.03 

Q go, +0.006 

U.U01_o 005 

8.8(f) 
6.97(f) 
056+0-019 

U.UOD_o 020 

n n'i9+o-oo5 

U.Ud9_o 002 

9.3(f) 
4370.1/3288 



Obs 2 



2 69+0-03 
^•"^-0.02 

9 35+0-1* 
"••-'''-0.06 

n «ifi+0-007 

0.»ib_o 006 

9 O+O-'^ 
"•'-'-0.6 
169+0 074 
u.iua_o 007 
p +0.002 
U.zoi_o 007 
9 7c;+0.03 
^•'^-0.02 



6.7(f) 

0.004 ± 0.001 

n r,r, + 0.08 

^•^^-0.07 

8.8(f) 
6.97(f) 
< 0.05 

013+0-023 
U.Uici_o 005 



9.3(f) 



Notes. -Model 1 is described in XSPEC as PHABSxKDBLURx(PL+REFHlDEN). The value of N-n, inclination and rin were tied between the two 
observations. An instrumental line at 1.876 keV was added to each model. The normalisation of each component is referred to as N. Frozen values are 
followed by (f). -ff^en i** given in units of 102l H cm~3 Model 2 is similar to the previous model and only fitted to the data below 6.6keV. Model 3 and 4 
includes Fe XXV and Fe XXVI absorptions (GABS model in XSPEC) at fixed energies of 6.7 and 6.97 keV respectively. Their respective absorption edges 
is modelled with the model EDGE in XSPEC with energies frozen at 8.8 and 9.3 kcV respectively. The optical depth r of the absorptions and edges are 

linked for consistency. 
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Figure 8. (a): vs plot for J1655-40. Using Model 4 (Table 2) with 
the inclination allowed to vary over its full range we obtain a value of 
2.17^0 j7'"g f™' the inner radius at the 90 per cent confidence level for one 
parameter of interest (A.x^ =2.71 criterion) shown by the solid horizontal 
line.(b): Similar plot for the inclination of J1655-40. 
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Figure 9. Best-fit spectra for J1655-40 using Model 4. Obs 1 and 2 are 
shown in black and red respectively. The complex absorptions are phe- 
nomenologically modelled with two Gaussian absorption (see text). The 
data have been rebinned for plotting puiposes only. 



values as low as 1.3rg (spin > 0.99) obtained with an inclination 
of 70 degrees (Model 2) rising to 2.32rg (spin > 0.9) with the in- 
clination allowed to vary (Model 4). However, it must be noted that 
using the reflection model both with and without accounting for 
the absorption features as well as both with and without constraints 
on the inclination, gives results with an inner-radius consistently 
below 2.32rg (Table 2). Taking this value as an upper limit to the 
innermost radius of emission implies a dimensionless spin param- 
eter greater than 0.9 for J1655-40 (Fig. 6). 
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Figure 10. Inclination-inner radius contour plot for J1655-40. The 68, 90 
and 95 per cent confidence range for two parameters of interest are sliown 
in black, red and green respectively. It can be seen that for the inner radius is 
constrained to be below approximately 2.5rgfor a large range of inclination 
at the 90 per cent confidence level for two parameters. 

4 DISCUSSION 

The X-ray spectra of galactic black hole binaries can provide im- 
portant information on both the geometry of the system as well 
as intrinsic physical parameters such as the spin of the central 
black hole. Using the reflection features embedded in the spec- 
tra of both J1753.5-0127 and J1655-40 we found clear indica- 
tions that the accretion disk extends close to the radius of marginal 
stability in both cases. For J1753. 5-0127 we have shown that 
the innermost emitting region extends down to rin= S.lto's'g 
(Fig. 4a) with an innermost inclination of SSly degrees (Fig. 
4b). Based on the normalisation of the disk blackbody compo- 
nent, Miller, Homan & Miniutti (2006) found an inner radius of 
nn= 2.O(3)(Af/lOM0)(d/8.5kpc)(cosi)"^/^rg, similar to ours 
for a large range of parameters. Our result is also consistent with 
that found by Soleri et al. (2008) of 2.6-6.0rg for a 10 M© black 
hole. In both these cases nn was estimated using the normalisation 
of the MCD component and thus requires knowledge of the mass, 
distance and inclination of the sources. Uncertainties in these val- 
ues can significantly affect the value of the inner radius and thus 
spin parameter. 

Estimating the spin of a black hole based on the extent of the 
emitting region relies on the assumption that the accretion disk ex- 
tends to the innermost stable circular orbit (ISCO) and that emis- 
sion within this radius is negligible. Reynolds & Fabian (2008) ad- 
dressed the robustness of this assumption and found that reflection 
within the ISCO becomes significantly less as one considers more 
rapidly rotating black holes. The way the dimensionless spin pa- 
rameter depends on the position of the ISCO is shown in Fig. 6 
(soUd curve). Assuming that the accretion disk extend to the ISCO 
for 11753.5-0127 we obtain a spin parameters of 0.76lo j5. Note 
that this estimate is independent of the unknown distance to the 
source. 

The system in J1655-40 is known to have an orbital-plane in- 
clination of approximately 70 degrees (Greene, Bailyn & Orosz 
2001). Constraining the inclination to this value resulted in an in- 
ner radius of emission consistent with that of a maximally rotating 
black hole (Model 3, Table 2). Allowing the inclination of the blur- 
ring function to be a free parameter in the spectral fit of 11655-40 
we obtain an improved fit with an inner radius of = 2.ntQ^^r^ 
(Fig. 8a) and an inner-disk inclination of 30^'iq degrees (Fig. 8b). 



This value for the inner radius implies a lower limit for the spin pa- 
rameter of 0.9 (Fig. 6a). The possibility of a misalignment between 
the innermost region of the accretion disk and the orbital plane 
inclination in galactic microquasars have been suggested by vari- 
ous authors (see e.g. Maccarone 2002) with clear examples seen 
in both 11655-40 (Martin, Tout & Pringle 2008) and V4641 Sgr 
(Martin, Reis & Pringle 2008). It was shown in Martin, Tout & 
Pringle (2008) that the alignment time scale in microquasar such as 
J 1655-40 is usually a significant fraction of the lifetime of the sys- 
tem. Thus if the black hole in such a system were formed with mis- 
aligned angular momentum, as expected from supernova-induced 
kicks, then it would be likely that the system would remain mis- 
aligned for most of its lifetime. This assumption is supported by 
the apparent misalignment found in the present work. 

The lower limit for the spin parameter of J 165 5 -40 found here 
(>0.9; Fig. 6) is not consistent with that reported by Shafee et al. 
(2006) of 0.65-0.75. As was mentioned above their method re- 
quires prior knowledge of several other factors including the inner 
disk inclination and distance to the source. Pszota & Cui (2007) 
have recently shown that neither disk continuum models KERRBB 
(Li et al. 2005) nor BHSPEC (Davis & Hubeny 2006) were able to 
successfully model the ultrasoft spectra of GX 339-4, nonetheless 
their best fit parameters in conjunction with the best estimates for 
the physical parameters of the source, suggested a moderate spin 
of 0.5-0.6 in comparison with our estimate of « 0.935 (Reis et al. 
2008; Miller et al. 2008) based on the various reflection features 
in the spectra of GX 339-4. Our suggestion that 11655-40 contains 
a rapidly rotating black hole is in agreement with results based on 
Quasi-Periodic Oscillations where a value of > 0.91 is usually de- 
rived (Zhang et al. 1997; Cui et al. 1998; Wagoner et al. 2001; 
RezzoUa et al. 2003). Furthermore, in their derivation of the spin 
parameter of J 1655-40, Shafee et al. (2006) assumed both a lack 
of misalignment between the central disk region and the orbital 
plane and more importantly they used a distance to the source of 
3.2 ± 0.2 kpc. Foellmi (2008) has recently shown that the distance 
to J1655-40 is most likely to be less than 2 kpc. When this distance 
is used in place of 3.2 kpc the spin parameter of J1655-40 derived 
by the thermal continuum method becomes greater than 0.91 in 
agreement with our results. 



5 CONCLUSIONS 

We have studied the XMM-Newton spectra of both SWIFT 
11753.5-0127 and GRO 11655-40. By estimating the innermost ra- 
dius of emission in these systems we constrain their spin param- 
eters by assuming that the accretion disk extends down to the ra- 
dius of marginal stability. For J1753. 5-0127 the spin is found to be 
0.761q ^5 at 90 per cent confidence. The innermost disk inclina- 
tion in J1753. 5-0127 is estimated at 55^7 degrees. In the case of 
11655-40 we find that the best fit requires a disk which is signifi- 
cantly misaligned to the orbital plane. An inclination of SOl^g de- 
grees and a dimensionless spin parameter greater than 0.9 is found 
at the 90 per cent confidence level. Our method involves spectral 
modelling of both the intrinsic disk thermal emission as well as 
the reprocessed hard radiation which manifest itself as various re- 
flection "signatures". These reflection signatures are independent 
of both black hole mass and distance and are thus a very useful tool 
for such measurements. 



8 R. C. Reis et al. 



ACKNOWLEDGEMENTS 

RCR acknowledges STFC for financial support. ACF and RRR 
thanks the Royal Society and the College of the Holy Cross re- 
spectively. 

Note added after submission: In a recent arXiv posting, Hiem- 
stra et al. (2009) published an analysis made on the same XMM- 
Newton and RXTE data for J1753. 5-0127 used here. They confirm 
the existence of a broad iron-A"a line and show that various contin- 
uum models with the addition of a LAOR line can successfully fit 
the data. More importantly, they claim to successfully fit the X-ray 
spectrum of J1753. 5-0127 with a continuum model that does not re- 
quire emission from a disk extending down to the ISCO. We have 
investigated this solution and find that in their case the broad line in 
a consequence of Compton broadening of the emitted photons by 
the hot surface of the accretion disk (Ross, Fabian & Young 1999; 
Ross & Fabian 2007). Their solution requires a highly ionised disk 
(^ ~ 5000 erg cm s^^) truncated at ~ 255 Tg with an inclina- 
tion angle consistent with zero. Using the same reflection model 
(REFLION; Ross & Fabian 2005) convolved with KDBLUR and 
including a disk component, we find a further solution where the 
disk extends down to the ISCO (6rg) and has a much lower ionisa- 
tion parameter of ^ ~ 500 erg cm s^^. This second interpretation 
yields an inclination of ~ 65 degrees, and an improvement of -7.7 
in for 2 extra degrees of freedom. The solution presented in our 
paper reflects that of the lower ionisation reported above. We are 
concerned that in a solution with such high ionisation parameter 
and a truncation radius at ^ 45 times the ISCO would result in 
a disk which is approximately 2 x 10'' times less dense that that 
of a similar disk extending down to the ISCO. With an estimated 
mass of approximately 10 M0 (Cadolle Bel et al. 2007), a radius of 
255rg is the equivalent of a disk truncated at ~ 3.8 x 10* cm from 
the central black hole. Using the ionisation parameter and unab- 
sorbed flux presented in Hiemstra et al. (2009) we can estimate the 
hydrogen number density to be approximately 2.3 x lO'^^cm^^. 
A modest surface layer of Thomson depth tt ~ 3 will thus re- 
sults in a disk with a half-thickness, t greater than 2 x 10^ cm; 
at least 5 times larger than the truncation radius. For a disk with 
inner radius r >>rg we can write the ionisation parameters as 
^ « 1.13 x 10^77 t/r^ergcms"^, where rj is the X-ray efficiency 
(Reynolds & Begelman 1997). Fig. 1 1 shows the ratio f/r as a func- 
tion of r for a system with 77 = 0.06 and an ionisation parameter 
of 5000 and 500 erg cm s^^. For a typical thin accretion disk t/r 
should be well below unity and as can be seen in Fig. 1 1 there are 
no solutions in this range for a disk with ^ ~ 5000 erg cm . We 
conclude that the high ionisation solution found by Hiemstra et al 
(2009) is physically inconsistent; the disc must extend in to small 
radii. 
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